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In a recent publication we established an analogy between the free energy of a hard sphere system 
and the energy of an elastic network 0|. This result enables one to study the free energy landscape 
of hard spheres, in particular to define normal modes. In this Letter we use these tools to analyze 
the activated transitions between meta-bassins, both in the aging regime deep in the glass phase and 
near the glass transition. We observe numerically that structural relaxation occurs mostly along 
a very small number of nearly-unstable extended modes. This number decays for denser packing 
and is significantly lowered as the system undergoes the glass transition. This observation supports 
that structural relaxation and marginal modes share common properties. In particular theoretical 
results 0,0] show that these modes extend at least on some length scale I* ~ (<f) c — 0)~ 1//2 where 
(f> c corresponds to the maximum packing fraction, i.e. the jamming transition. This prediction is 
consistent with very recent numerical observations of sheared systems near the jamming threshold 
jfjj, where a similar exponent is found, and with the commonly observed growth of the rearranging 
regions with compression near the glass transition. 

PACS numbers: 

A colossal effort has been made to characterize the spatial nature of the structural relaxation near the glass 
transition. Numerical simulations [B[ and experiments @, Q have shown that the dynamics in super-cooled liquids 
is heterogeneous. Both the string-like [|| and the compact 0] aspects of the particles' displacements have been 
emphasized. Nevertheless, the cause of such collective motions remains debated [13, El- To make progress, one would 
like to relate these motions to other objects. A possible candidate is the excess of low-frequency modes present in all 



glasses, the so-called boson peak (12(. Because these modes shift in general to lower frequencies as the temperature 
increases toward the glass transition temperature T g , it has been proposed that they are responsible for the melting 



of the glass [13|, UM- This suggests the use of widely employed tools, such as the low- frequency instantaneous normal 
modes 15 or the negative directions of saddles of the potential energy landscape [16[, to analyze the collective 
motions causing relaxation. Nevertheless, this approach has the major drawback of being based on energy instead 
of free energy. As such, it cannot be applied for example to hard spheres or colloids, where structural relaxation is 
also known to be collective, see e.g. [7|. In this case barriers between meta-stable states are purely entropic. More 
generally, one expects entropic effects to be important for glasses where hard-core repulsions and non-linearities are 
not negligible, which is presumably the case in general above T g fl7l ] . 

Recent developments make this analysis possible in hard sphere systems. In [l[, we derived an analogy between the 
free energy of a hard sphere glass and the energy of a weakly-connected network of logarithmic springs. This allows 
us to define normal modes, that can be compared with the dynamics. Furthermore, recent results 0, H [13 valid 
for weakly-connected networks, such as clastic particles near jamming (lil. l2Cl] — where scaling laws between packing 



geometry and vibrational properties were first observed — or simple models of silica [21( , apply to characterize these 
modes: (i) Excess modes appear above some frequency w AM which depends on the pressure p and the coordination z, 
whose definition shall be recalled below for hard spheres. These anomalous modes extend at least on a length scale 
/*, which depends on z and diverges near maximum packing 

BE, El- ( u ) Meta-stable states can exist only if they 
contain a configuration for which uj am > 0. This leads to a non-trivial scaling relationship between p and z that must 
be satisfied in the glass phase. Numerically, we observed that the hard sphere glass lies close to marginal stability: 
the coordination is just sufficient to maintain rigidity [![. This implies that anomalous modes are present at very 
low-frequency 

In this Letter, we study how low-frequency modes take part in the structural relaxation, both during the aging 
dynamics deep in the glass phase, and in the vicinity of the glass transition where the system is at equilibrium. We 
show that when relaxation occurs between meta-stable states, the system yields in the direction of the softest modes: 
most of the amplitude of the observed displacements can be decomposed on a small fraction of the modes, of the 
order of few percent. This observation supports that the collective aspect of the relaxation does not stem from the 
non-linear coupling of localized relaxation events, but rather from the extended character of the softest degrees of 
freedom. This suggests that the typical size of the events relaxing the structure increases as the extension of the 
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anomalous modes I* ~ (cj> c — 0) -1 ^ 2 ~ p 1 ^ 2 , which diverges deep in the glass phase. 

We start by recalling some results of [l[ . In a meta-stable state of a hard sphere system, one can define a contact 
network [23[ : two particles are said to be in contact if they collide during some interval of time t\ , where ti is chosen to 
be much larger than an r c , the collision time, and smaller than the structural relaxation time r where meta-stability 
is lost. The coordination number z of this network is defined as the average number of contacts of the particles in the 
system. An approximation of the Gibbs free energy Q can then be expressed as a sum on all the contacts (ij): 

g = -wrYMMt) (i) 

(ij) 

where hij = m — r- t — rj is the gap between particles i and j, rij is the distance between them, denotes the radius of 
particle i, and ()t is a time-average. Eq.([T]) has two main limitations: (i) it is only exact near the maximum packing 
fraction <j> c where the pressure diverges and (ii) to perform the time-average one requires a strong separation of time 
scales between r c and r. Thus Eq.© is a better approximation deep in the glass phase. Nevertheless the corrections 
to Eq.© are found to be rather small empirically [H, HH, and we shall use Eq.JT]) to study the vicinity of the glass 
transition {4> m <J)q) as well. 

Eq.© can be expanded around any equilibrium position. For a contact ij, one finds for the force = — kT / (hij)t 
and for the stiffness Vy = kT/(hij)^. This enables one to compute the dynamical matrix M. [27J which relates a 
small applied force to the linear displacement of the average particle positions. Normal modes can then be computed, 
whose angular frequencies are the square roots of the eigenvalues of Ai. In what follows we locate quiet periods 
of the dynamics where Ai can be estimated. Then, we use the normal modes to analyze the subsequent structural 
relaxation. 

We consider a bidisperse two-dimensional hard sphere system. Half of the particles have a diameter o\ = 1, the 
other a diameter 1.4, their mass is m = 1, and energies are expressed in units of kT. To study the aging dynamics, 
configurations are generated in the glass phase (</>o ~ 0.79 < < <p c ~ 0.84) as in [l[. An event-driven code is used to 
simulate the dynamics. We observe long quiet periods, or meta-stable states, interrupted by sudden rearrangements, 
or "earthquakes" . Such earthquakes correspond to collective motions of a large number of particles, and have been 
observed in various other aging systems, such as colloidal paste or laponite [261 ] . and in Lennard- Jones simulations 
[2^. [25I] . Even for our largest numerical box of N = 1024 particles, deep in the glass phase these events generally span 
the entire system. They appear as drops in the self scattering function C(q,t) = (exp[iq.(Ri(t) — i?j(0))])j, where ()j 
is an average on all particles and Ri(t) is the position of particle i at time t. An example of earthquake is shown in 
Fig.©. 

In what follows the average particles position in a meta-stable state / is denoted \R l ) = {(Ri)t},i = 1...N. In 
practice the time-average ()t is done on a long time t\ corresponding to few hundred collisions per particle (we use 
t\ = 10 5 , 5 x 10 4 numerical time steps for respectively N = 1024 and 256 particles). The earthquake displacement 
field \SR e ) between two meta-stable states I and m is then defined as \SR e ) = \R m ) — \R l ), see Fig©b). During 
earthquakes, we find that the average particles displacement is typically 10%-20% of the particles diameter, and tend 
to decrease with the pressure. 

To analyze these displacement fields, we compute the average of the particles' positions and the contact network in 
the meta-stable state prior to the earthquake[33|. This enables us to define M. and the normal modes \SR a ) : where 
the label a = 1, ...,2N ranks the modes by increasing frequencies w° . An example of the density of states D(u>) 
is shown in Fig.([2]-a), together with the lowest frequency mode, Fig.([2]-b). We indeed observe extended anomalous 
modes at very low frequencies, in agreement with the marginal stability inferred from the microscopic structure of 
the glass [l[. Note that we occasionally observe a few unstable modes even deep in the glass phase, implying the 
presence of saddles (and multiple configurations of free energy minima) or "shoulders" in the meta-stable states we 
are considering. In the present work we do not focus on this aspect, and treat unstable modes as the rest. 

We then project the earthquake displacement \SR e ) on the modes and compute c a = {SR e \SR a ) / (5R e \SR e ) , where 
{5R e \5R a ) = J2i ' ■ The c Q 's satisfy J2a °a 2 = 1 since the normal modes form a unitary basis. To study how 
the contribution of the modes depends on frequency, we define g(u>) = (c Q 2 ) w , where the average is made on all a 
such that u! a e [uj,lu + duj}. Fig.([2]-a) shows g(ui) for the earthquake shown in Fig.©. The average contribution of the 
modes decreases very rapidly with increasing frequency, and most of the displacement projects on the excess-modes 
present near zero-frequency. This implies that the free energy barrier crossed by the system is located in the direction 
of the softest degrees of freedom. 

To make this observation systematic, we introduce the label i to rank the c's by decreasing order: c\ > C2... > C2N- 
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FIG. 1: Left: Self-density correlation function C(q,t) vs. time for q — 2tt/<j 1 in a system of N — 256 particles, at packing 
fraction cj> = 0.837. Meta-stable states appear as plateaus of C(q, t), whereas the drops of C(q, t) are the mentioned earthquakes. 
Time-averages are made during the time segments t\. Right: Displacement field of the corresponding earthquake. Arrows 
connect the average particles position before and after the earthquake, they are amplified 4 times here for visibility. For similar 
data in a 3d LJ see 1241 1. 




FIG. 2: Left, straight curve: D(uj) computed in the meta-stable state prior to the earthquake of Fig.JTJ) vs. ui/(f), the angular 
frequency rescaled by the average contact force {/), computed by averaging on all pairs of particles in contact. Dotted curve: 
g(uj) vs. cu/{f). Right: Lowest frequency mode. 



Then we define: 

k 

F(k)^Y, c * ( 2 ) 



i=i 



F(k) indicates which fraction of the total displacement is contained in the k most contributing modes. If F(k) = 1 Vfc 
then only one mode contributes. If F(k) = k/2N all modes contribute equally. We then define ki/ 2 as the minimum 
number of modes contributing to 50% of the displacements, i.e. the smallest k for which F(k) > 1/2. Fig.([3|) shows 
F(k) and F 1/2 = k 1/2 /(2N) for the 17 cracks studied. Fig.([3}b) shows that 0.2% < F 1/2 < 2% for all the events 
studied throughout the glass phase. We thus systematically observe that the extended earthquakes correspond to the 
relaxation of a small number of degrees of freedom. 

We now extend this analysis to an equilibrated super-cooled liquid. We equilibrate at 0.77 < <f> < 0.786. As 
previously observed e.g. in [9j], the dynamics is heterogeneous in time, and sudden rearrangements still occur on time 



scales of the order of r, the time scale of the a-relaxation [34( . We use the procedure previously described to determine 



the average configuration of meta-stable states, and to define the displacement relaxing the structure, see Fig.((3]-a). 
We start by studying five packing fractions in a system of N = 64 particles. For this size rearrangements generally 
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FIG. 3: Left: Examples of F(k) vs. k/2N for systems with different average contact force {/). (/) is measured before the 
earthquake. It is proportional to the pressure near <f> c , and is of the order of 20 near the glass transition (f>o- Right: -Fi/2 vs. 
(/) for N = 256 (circles) and N — 1024 (diamonds) particles. 



span the entire system. A similar observation was made in a 3-dimensional lennard- Jones system of 125 particles [9(. 
For each packing fraction, F(k) and -F\/2 are computed for seven relaxation events. Results for (F 1 / 2 ) are presented 
in Fig.^b) as a function of packing fraction. We find that (F 1 / 2 ) < 4% for all <fi studied, implying that also in 
this region of the super-cooled liquid phase a small fraction of the low-frequency modes contribute to the relaxation 
events. Interestingly, this fraction decays significantly as the relaxation time grows, suggesting a rarefaction of the 
number of directions along which the system can yield near the glass transition. For the largest (f> « 0.786 we studied, 
(F1/2) < 2%, which implies that the collective event relaxing the system corresponds mainly to one or two modes. 
Fig.([5]-a) and Fig.((5]-b) compare one event and the mode that contributes most to it, which turns out to be the softest 
mode in this particular example. 

Size effect: We now consider a system of TV = 256 particles. An example of relaxation is shown in Fig.([5]-c). We 
observe for this system size that a larger fraction of the particles stay motionless. It is also interesting to compare for 
the same system size the aging dynamics deeper in the glass phase, e.g. Fig.^b) with Fig.([5]-c): as was previously 
observed in LJ systems 24j, the collective rearrangements at equilibrium involve less particles than earthquakes, but 
move them more. Nevertheless in the equilibrated case as well, we shall see now that the observed displacements 
projects on a very small fraction of the normal modes. To study this question we perform the analysis introduced 
above, and compute (-Fi/2} by averaging on twelve events for the five packing fractions considered. Results are 
shown in Fig.([4jb): they are qualitatively similar to the system with N = 64, but the values of {F 1 / 2 } are larger 
by approximately0.5% . Nevertheless, it is well know that the glass transition occurs at lower packing fraction for 
smaller systems [29(, see inset of Fig.(j4]-c). When plotted versus relaxation time, the difference between the values of 
(F1/2) m the two systems is reduced and less systematic, and (-F1/2} becomes in fact smaller for N = 256. Thus, even 
when embedded in a system containing quiet regions, relaxation occurs along the softest modes. 

In kinctically constrained models [30j, that have been proposed as paradigms for glassy dynamics, heterogeneous 
dynamics can arise from simple local microscopic rules. Our result that collective rearrangements correspond mostly 
to a few (and therefore necessarily reasonably extended) modes supports an alternative view: elementary relaxations 
events are already extended objects, as are the soft degrees of freedom of the system. We now argue that these 
modes are the anomalous modes 0, 9 described in the introduction. Our justification for this lies in the microscopic 
structure of the glass: as evoked above, imposing the marginal stability of these modes leads to a power-law relation 
between coordination and pressure that is indeed observed in the glass phase If other soft objects (e.g. local 
configurations particularly soft due to disorder) were driving the dynamics, it is unlikely that the glass would freeze 
in this specific region of the phase space [35| . 

More needs to be known on the statistical properties of the anomalous modes, for example on the curly aspect of 
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FIG. 4: (a): C(q,t) in an equilibrated system of iV — 64 particles at (f> = 0.786. The segment ti corresponds to 10 4 numerical 
time steps, which corresponds roughly to 300 collisions on average per particle, (b): {-F1/2) vs (f> for t\ = 5 X 10 4 (circle) for 
A*' = 256 and t\ = 10 4 (triangle) for a system with ./V = 64. The average number of collisions per particle is written in the 
legend and denoted by nc. (c): (-F1/2) vs t for two different system size. Inset: relaxation time r vs <j> for JV = 64 (triangle) 
and N = 256 (circle). 
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FIG. 5: (a) Displacement field corresponding to the relaxation of Fig.Q-a). Arrows were multiplied by 1.2. (b) lowest frequency 
normal mode, which has the highest projection in this particular case. In this example one mode contributes to most of the 
displacement, (c) Relaxation event in system of N — 256 particles at <f> = 0.786. (d) Projection of the displacement (c) on the 
3% of the mode that contribute the most to it. 



their displacement field or on their apparent capacity to form intense flow lines or strings. This problem turns out 
to be equivalent to the statistics of force chains in amorphous solids [3l1 |. a much studied problem in the granular 
literature, but which still lacks a definite answer. Nevertheless, we know that the anomalous modes are characterized 
by some length scale I* ~ p 1 / 2 above which the softest modes must extend 0. The observation that the softest modes 
dominate the relaxation supports that regions of size at least I* ~ p 1 ^ 2 ~ (<j> c — 0)" 1 / 2 must rearrange. Very recent 
numerical work of sheared system near them jamming threshold Q supports our views, as in this case rearrangements 
are characterized by a diverging length scale with an exponent 0.6+/- 0.1, in agreement with our predictions. These 
are also consistent with the growing dynamical length scale observed near the glass transition. Nevertheless in this 
case the length scales at play, typically about 5 or 10 particles, are too small to compare different theories (32|. More 
stringent tests could be performed near maximum packing, where a diverging length scale is expected. This could be 
tested experimentally, e.g. by considering the intermittent aging dynamics of colloids at large osmotic pressures. 

Our work supports a unified description of both the structural relaxation and the packing geometry, where the 
dynamics corresponds to the collapse of anomalous modes, and where the microscopic structure is fixed by their 
marginal stability. Note that the theoretical framework used here applies identically in three dimensions, where we 
expect our results to be valid as well. This scenario may also hold in other glasses, for example in Lennard- Jones 
where anomalous modes are also present 22 1. Nevertheless in this case, as for any long-range interaction potentials, 
I* is bounded above and does not diverge in the glass phase. 
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